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Abstract. General relativity has predicted the existence of gravitational waves (GW), which are 
waves of the distortions of space-time with two degrees of polarization and the propagation speed 
of light. Alternative theories predict more polarizations, up to a maximum of six, and possible 
deviation of propagation speed from the light speed. The present paper reviews recent proposals to 
test the gravity theories in the radiation regime by observing GWs using pulsar timing arrays. 
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INTRODUCTION 

Two important characteristics of GWs are important to differentiate the validity of grav- 
ity theories in the radiative regime; the polarization and dispersion of GW in vacuum. In 
alternative metric theories, GW can have up to six possible polarization states, four more 
then are allowed by GR, Furthermore, the propagation speed of GW can deviate from 
the predication of GR that GW propagates at light speed in vacuum, i.e. the effective 
graviton mass is zero. 

Pulsar timing array is a unique technique to detect nano-Hertz GWs by timing mul- 
tiple millisecond pulsars, which are very stable celestial clocks [1]. It turns out that 
a stochastic GW background leaves an angular-dependent correlation in pulsar timing 
residuals for widely spaced pulsars [2]: the correlation C{6) between timing residual of 
pulsar pairs is a function of angular separation 6 between the pulsars. One can analyse 
the timing residual and measure such a correlation to detect GWs [1]. Lee et al. (2008, 
2010) have found that the exact form of C{6) is very different from the one of GR, if 
the GW has extra polarization state or graviton mass is not zero. Thus by measuring the 
correlation function, we can directly test gravity theories in the radiative regime. 



PULSAR TIMING CORRELATION FUNCTIONS AND GRAVITY 

TESTS 



A GW introduces extra signal in pulsar timing data. Let the unit vector of the GW 
propagation direction be e^, GW frequency be /, the direction from the observer to the 
photon source (pulsar) be n'. The GW induced frequency-shift of a pulsar timing signal 



is [3, 4, 5] 
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where the D is the displacement vector from the observer to the pulsar, hij{t^O) and 
hij{t — |D|/c,D) are the metric perturbations by GW at the Earth and at the pulsar 
when the received pulse was emitted, (O is the angular frequency for the pulsar pulse, 
/cut = rugC^ /h is the cut-off frequency of GW due to the graviton mass mg, and the is 
the GW wave vector given by [4] 

kg(/)= ^^^^';^'"^^% .. (2) 

The induced pulsar timing residuals R{t) are given by the temporal integration of above 

the frequency shift at Earth given above, thus R{t) = Jq ^^^dT. 

The spatial metric perturbation hij{t,r) induced by a stochastic GW background is a 
superposition of monochromatic GWs with random phases and amplitudes. It is [6] 

hij{ty)= E r df [ dah^'ifAWji^y^^''^'-^^^^^^ (3) 
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where Q. is the solid angle, index /, j run from 1 to 3, is the amplitude of the GW 
propagating in the direction of per unit solid angle per frequency of polarization state 
P, and the polarization tensor e^^ of GWs are given in details in Lee et al. (2008). The 
superscript P takes value of '+, x' for the two Einsteinian modes of GW polarization, 
and 7' for the breathing and the longitudinal mode respectively, and '5/1, se' for the 
shear modes. 

Such stochastic GW background leaves a correlation between timing residuals of 
pulsars pairs [2, 3]. Such correlation, C(0), depends on the angular distance Q between 
two pulsars as well as on the polarization of GW and graviton mass. 

Lee et al. (2008) have calculated the pulsar timing correlation function for all the 
polarization modes of GW. For the Einsteinian modes and for the breathing mode, the 
cross-correlation function C^(0) is independent of earth-pulsar distances and indepen- 
dent of the GW characteristic strain spectrum. In contrast, for the modes that are not 
purely transverse, the shear and longitudinal modes, the cross correlation functions de- 
pend on the specifics of the strain spectra and on the pulsar distribution in distance. 

Fig. 1 shows the correlation function according to different classes of GW polariza- 
tion. Clearly by comparison of these 'theoretical' correlation curves with observations 
we can test the polarization state of GWs. 

Lee et al. (2010) have calculated the pulsar timing correlation function for a GW 
background with none-zero-mass graviton. They noted that the pulsar timing cross- 
correlation function for a massive GW background depends on the graviton mass, 
specific power spectra of the GW background, and on the observation schedule. The 
5-year and 10-year correlation functions are reproduced in Fig. 2, where the graviton 




FIGURE 1. The normalized pulsar timing residual correlation coefficients. Here 9 is the angular 
separation between two pulsars. 'GR' stands for the two transverse traceless modes, ' +' and ' x ' . Results 
are given for several values of a, the power-law index of the GW spectrum. 
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FIGURE 2. The atlas for cross-correlation functions C{6). The label of each curve indicates the 
corresponding graviton mass in unit of electron-volts (eV). The left panel are the correlation functions 
for a 5-year bi-weekly observations. The right panel shows correlation functions for 10-year bi-weekly 
cases. We take a — —2/3 for these results. These correlation are normalized such that the C(0) = 0.5 for 
two different pulsars. 
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with the same mass introduces more deviation to the 10-year correlation function than it 
does to the 5 -year one. 

Intuitively speaking, the necessary conditions for a positive detection of a graviton 
mass should be: 1. The GW is strong enough such that the GW can be detected; 2. the 
physical effects of altemative theories should be strong enough to see the deviation from 
GR. These intuition is confirmed by simulations in [3, 4], which show that the high 
detection rate is achieved only if one has enough pulsar and if the graviton mass is large 



enough or if GW of alternative polarization modes is strong enough. 

For identifying the polarization modes, observation shows that if bi-weekly observa- 
tions are made for five years with RMS timing accuracy of 100 ns, then 60 pulsars are 
required for the longitudinal mode; 60 for the two spin-1 "shear" modes; and 40 for the 
spin "breathing" mode and 40 pulsars are needed for the detection of the GR modes. 

For detecting massive graviton, simulations have shown that we need at least 60 
pulsars to be able to tell the difference between a massive GW background and a 
massless one. For 5-year timing of 100 pulsar we can start to detect a graviton heavier 
than 2.5 x 10^^^ eV and we can achieve a limit of nig = 10^^^ eV by using 5-year 
observation of 300 pulsars. We can achieve levels of 10^^^ eV and 5 x 10~^^ eV in a 
10-year observation using 100 and 300 pulsars respectively. 

CONCLUSION AND DISCUSSION 

The stochastic GW background produces extra timing signals in pulsar TOA data, and 
one can detect the GW background by precise timing several pulsars and measuring the 
angular dependent correlations between the timing signals of several pulsars. A precise 
measurement of the angular correlation function can, in principle, determine the GW 
polarization properties of the GWs making up the stochastic background as well as the 
graviton mass. A large number of pulsars with good timing precise of ~ 100 ns are 
required to successfully perform tests for gravity theories in radiative regime. In this 
regard, pulsar surveys' success in finding more millisecond pulsars is critical. To time 
such large mount of pulsars, the Large European Array for Pulsars [7], the Five-hundred- 
meter Aperture Spherical Radio Telescope [8, 9] and the Square Kilometer Array (SKA) 
will offer unique opportunities to detect the GW background and measure its properties. 
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